Jahn- Teller Distortion in Bimetallic Oxalates 
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A Cs-symmetric crystal- field potential in the Fe(II)Fe(III) bimetallic oxalates splits the L — 2 
Fe(II) multiplet into two doublets and a singlet. In compounds that exhibit magnetic compensation, 
one of the doublets was predicted to lie lowest in energy and carry a non-quenched orbital angular 
momentum ±L^-^, where exceeds a threshold value. In a range of , a Jahn- Teller (JT) 
distortion increases the energy splitting of the low-lying doublet and breaks the C3 symmetry of 
the bimetallic planes around the ferrimagnetic transition temperature. At low temperatures, the 
JT distortion disappears in compounds that display magnetic compensation due to the competition 
with the spin-orbit coupling. A comparison with recent measurements provides strong evidence for 
this re-entrant, low-temperature JT transition and a prediction for the normal, high-temperature 
JT transition. The size of the JT distortion is estimated using first-principles calculations, which 
suggest that the long-range ordering of smaller, non-Cs-symmetric organic cations can eliminate 
magnetic compensation. 



I. INTRODUCTION 

Jahn- Teller (JT) transitions [1] in which electronic de- 
generacies are removed by crystal distortions have been 
observed quite frequently in recent studies of molecule- 
based magnets [2]. It is well-known that the JT transi- 
tion will be quenched if the degenerate levels carry or- 
bital angular momentum and the spin-orbit coupling is 
sufficiently strong. When the spin-orbit coupling and JT 
energies are comparable, however, a pseudo-JT transi- 
tion with rather interesting behavior is possible. In this 
paper, we study the JT transition in a class of molecule- 
based magnets where the spin-orbit coupling can be mod- 
ified by choosing different organic cations to lie between 
bimetallic layers. Due to the competition between the 
spin-orbit and JT energies, the JT distortion may vanish 
at low temperatures in a re-entrant, first-order transition. 

One of the most fascinating classes of molecule-based 
magnets, bimetallic oxalates A[M(II)M'(III)(ox)3] were 
first synthesized [3] in 1992. In the open honeycomb 
structure of each bimetallic layer, sketched in Fig. 1(a), 
the transition- metal ions M(II) and M'(III) are coupled 
by the oxalate bridges ox=C204 [4]. Depending on the 
metal atoms, a single bimetallic layer can be either fer- 
romagnetic or ferrimagnetic (M(II) and M'(III) moments 
parallel or anti-parallel) with magnetic moments pointing 
out of the plane. While the organic cation A separating 
the layers cannot alter the sign of the exchange coupling, 
it does affect the overall properties of the system. 

Some Fe(II)Fe(III) bimetallic oxalates exhibit mag- 
netic compensation (MC) due to cancellation of the mo- 
ments on the Fe(II) and Fe(III) sublattices below the fer- 
rimagnetic transition temperature [5] . Based on sym- 
metry and energy considerations, we recently explained 
why MC occurs for some organic cations but not for oth- 
ers [6]. The Ca-symmetric crystal field (invariant un- 
der in-plane rotations of 27r/3) splits the L = 2 Fe(II) 
multiplet into two doublets and a singlet. By shifting 
the Fe(II) ions with respect to the oxalate molecules, 
the cation A determines the average orbital angular mo- 



mentum ii^^*^ carried by the low-energy doublet. Com- 
pounds exhibit MC when the low-energy doublet lies be- 
low the singlet and L^^ exceeds the threshold Ic- For com- 
pounds that do not exhibit MC ("normal" compounds), 
either L^^ < Ic or the singlet state lies lowest in energy. 

In the absence of spin-orbit coupling or when L^J = 0, 
a JT distortion that breaks C3 symmetry will always 
split the low-energy doublet. For large enough L^-^, the 
spin-orbit coupling will quench the JT distortion. How- 
ever, in a range of L^J that includes /c, we obtain a JT- 
distorted phase between temperatures T^r^ and that 
bracket T^. Both the re-entrant, low-temperature tran- 
sition Tf^ and the normal, high-temperature transition 

are first order. Comparison with recent experiments 
[7] allows us to estimate the normal JT transition tem- 
perature in MC compounds. 

This paper is divided into 5 sections. Section II devel- 
ops the phenomenological model for the JT transition, 
with results described in Section III. First-principles cal- 
culations are described in Section IV and a discussion is 
contained in Section V. 



II. MODEL FOR THE JAHN-TELLER 
TRANSITION 

Our model assumes a hierarchy of three energy scales. 
The dominant energy is the Hund's coupling that de- 
termines the spins 5* = 2 and S' = 5/2 on the Fe(II) 
(3d^) and Fe(III) (3d^) sites. Next in importance is the 
Ca-symmetric crystal field V generated by the 6 oxygen 
atoms around each of the Fe sites. These 6 oxygen atoms 
form two triangles of slightly different sizes rotated by 
about 48 degrees with respect to each other, one above 
and the other below the plane of Fe atoms. The small- 
est energies are the antiferromagnetic exchange coupling 
JcS • between the Fe(II) and Fe(III) moments within 
each bimetallic layer, the spin-orbit coupling AL • S on 
the Fe(II) sites (A ^ -12.65 meV [8]), and any non-Ca- 
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FIG. 1: (Color online) (a) A portion of the open honey- 
comb lattice, displaying three equivalent displacements of the 
Fe(II) ion into the adjacent hexagons, (b) Two possible level- 
splitting schemes of the 3d^ Fe(II) orbitals in a Cs-symmetric 
crystal- field, with either the doublet -01, 2 or the singlet 1^3 
lowest in energy. 



symmetric contributions of the crystal potential. 

As shown in Fig. 1(b), the Ca-symmetric potential V 
splits the L = 2 multiplet of the Fe(II) sites into two 
doublets ipi^2 and '^4,5, and a singlet ips- The low-energy 
doublet V^i,2 carries an average orbital angular momen- 
tum -^L^J ^ where L^J ranges from to 2 depending on 
the crystal-field parameters. Of course, the singlet state 
'03 does not carry any orbital angular momentum. 

Including the orbital contribution, the Fe(II) magnetic 
moment can be written M(T) Ms{T) + Ml{T) < 0, 
where Ms = 2{Sz) and Ml = {lJ). Of course, the 
Fe(III) magnetic moment M'{T) = Ms'{T) = 2{S'^) > 
has no orbital contribution. The mean-field (MF) ap- 
proximation is used to treat the nearest-neighbor ex- 
change 

JcS • ^ Jc{ + S,Ms^)/2 - MsMs^/4} (1) 

between neighboring Fe(II) and F(III) spins. The MF 
eigenvalues of ipi,2;a and Tp^^ are 



ei, = i-\X\L'=/ + 3JcMs'/2)a, 



e2<T 



ilML'^/ + 3JeMs'/2)a, 



= {3JeMs'/2)a + A, 



(2) 



(3) 



(4) 



SO that the doublet '0i,2;cr is split by the spin-orbit cou- 
pling with |e2cr — eicr| = 2\Xa\Ll^ . Because the crystal- 
field potential is assumed much larger than the spin-orbit 
and exchange energies, |A| Jc or |A|. 

Using Eqs.(2) and (3) for the eigenvalues €1^2- a ^ it is 
straightforward to evaluate M{T), M'{T)^ and the av- 
erage magnetization M''^^{T) = {M{T) +M'(T))/2 = 



{\M'{T)\ - \M{T)\)/2. When 1^ ^ 0.238 < L^^^ < 1, 
M^^^iT) passes through at the compensation temper- 
ature Tcomp < Tc. Just below Tc, M^^^{T) < be- 
cause the parallel spin-orbit coupling between L and S 
causes the Fe(II) moment M{T) to increase in amplitude 
more rapidly than the Fe(III) moment M'{T). But at 
zero temperature, the Fe(III) moment saturates at the 
larger value |M'(0)| = 2S' > |M(0)| =25' + L""/ so that 
M^^s(O) > 0. 

After comparing the observed values of Tcomp /^c ^ 
0.62 and Tc ~ 45 K with the theoretical predictions, we 
estimated [6] that Jc ^ 0.46 meV and L""/ ^ 0.28 in 
MC compounds. Normal compounds can fall into two 
categories. Either the doublet remains lower in energy 
than the singlet (A > 0) but with L""/ < l'^ ^ 0.234 or 
the singlet lies lowest in energy (A < 0). Our model also 
predicted that two compensation points were possible in 
the narrow window l'^ < L^^ < 1^ The recent observation 
by Tang et al. [7] of two compensation points in the 
compound N(n-C4H9)4[Fe(II)Fe(III)ox3] (data shown in 
Fig. 3) would seem to confirm this prediction. 

If a doubly-degenerate level is occupied by a single elec- 
tron in the absence of spin-orbit coupling, then a local 
JT displacement the Fe(II) ions corresponding to one of 
three equivalent directions is always favored by the elec- 
tronic energy. The eigenstates ?/^icr and ?/^2cr on the Fe(II) 
sites are mixed by the JT distortion as described by the 
Hamiltonian [9] 



(5) 



where ^ is independent of a. If eicr = e2cr are given by the 
MF result eo^ = 3JcMV/2, then the eigenstates oi H^''' 
are ^jjaa = (l/V2)(V^ia + i^2a) and ipb^ = (1/^2) (V^i^ - 
V^2a), with La,b;a = {i^a,b;a\L z\il^a,b;a) = and eigenvalues 
ea,6;cr = eocr ± So in the absence of spin-orbit coupling, 
the orbital angular momentum is quenched by the JT 
distortion and the doublet tl^a,b;a is split by 2|^|. The 
T = energy is then given by ^E/N = -SJcSS' - ^ + 
where the second term is an elastic restoring 
potential and N is the number of Fe(II) or Fe(III) sites 
per bimetallic layer. Since a > 0, the T = equilibrium 
value for ^ is |A|/2a > 0. 

Including the spin-orbit interaction AL • S, the eigen- 



values of H^^^ are 



eOcT + to- and Cba = ^Oa — to 



where t^j = — sgn(cr) y {\L%^ + Hence, the doublet 
splitting 5 = 2\tcr\ is enhanced by the JT effect (strictly 
speaking, the pseudo-JT effect when eicr 7^ e2cr)- Since 
there is no spin-orbit coupling when <j = 0, to = ^ and 
Lao = T50 = 0. For cr 7^ 0, the spin-orbit coupling main- 
tains a nonzero orbital angular momentum Laa = —Lba 
even in the presence of a JT distortion: 



so that L„„ = La 



XL"/ a + to 



{XL^J af + XLV + ^ 



(6) 



When ^ 7^ 0, the JT distortion 
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suppresses the absolute values \Lacr\ < for a = ±1 
and ±2. 

With a correction to avoid double counting, the MF 
free energy can be written 



-Tlog{z„Z„ie3^^^-^-'/^^} 



F _ 



where 



Z„ = 2 ^ g-3J.M,,a/2T cosh(t,/T), 



^iii = 2 ^ ( 



(7) 



(8) 



(9) 



The sums in the partition functions Z\\ and Z\\\ are over 
cr = 0, ±1, ±2 and a' = ±1/2, ±3/2, ±5/2. The second 
line in Eq.(7) corresponds to the elastic energy. 

Breaking C3 symmetry, a local JT distortion with ±^ 
involves the displacement of an Fe(II) atom either into 
one of the three open hexagons or towards one of the 
three neighboring Fe(III) atoms. The former distortions 
are sketched in Fig. 1(a). The anharmonic 73(^/| A|)^ term 
in Eq.(7) reflects the different energy costs for those two 
types of distortions. So the anharmonicity and source of 
the first-order JT transition arise quite naturally on the 
open honeycomb lattice. The sign of 73 does not affect 
any physical results and is chosen to be negative only 
so that C > 0. Any fluctuations between the distorted 
atomic configurations are assumed slow compared to the 
electronic time scales. 

It is simple to obtain the equilibrium values for Ms 
and Ms' from the extremal conditions dF/dMs = 
dF/dMs' = 0. Because the spin-orbit energy AL • S 
is treated exactly. Ml is not a variational parameter in 
the MF free energy and must be determined separately 
from the condition 
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Ml 



j- La,e-^'^^s'<^/2T sinh(t JT). (10) 



Of course, the equilibrium value for the JT parameter ^ 
is obtained by minimizing F with respect to ^. 



III. MODEL RESULTS 



FIG. 2: (Color online) The (a) average magnetization and (b) 
the JT mixing energy ^ normalized by |A| versus temperature 
T/|A| for a range of L^^ in increments of 0.01 using Jc/\\\ — 
0.037 and the elastic constants given in the text. 



The average magnetization and JT energy ^ are plot- 
ted versus temperature in Figs. 2(a) and (b) for a — 3.7, 
73 = —1.9, and 74 = 1.1. As expected, the MC threshold 
Ic ^ 0.282 is enhanced by the JT distortion. When L""/ < 
0.260, the JT distortion ^(T) persists down to T = 0. For 
L^J^ = (but stih taking A > 0), ^(0) ^ 0.842|A| ^ 10.5 

meV and rj^^ ^ 0.578|A|. For L""/ = 0.2, the JT dis- 



tortion would increase the splitting of the a 
from 4|A|L^'^ = 10 meV to 2x1 {2X1%^ f 



= 2 doublet 
^2 ^ 23 meV. 



When L^J > 0.260, the JT distortion is quenched at 
T = due to the strong orbital ordering. Hence, we ob- 
tain both lower and upper JT transitions, Tj!^ and Tj!^ . 
At the re-entrant JT transition, the distortion ^ vanishes 
below Tj!^ and appears above t|x- "^^^ temperature 



Tj!^ decreases as L^J increases and vanishes 



range ijrp^ 

when L""/ > 0.324. Notice that the JT transitions at t]^^ 
and Tj!^ always bracket Tc. 

Because the orbital contribution Ml{T) to the Fe(II) 
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FIG. 3: (Color online) The average magnetization versus 
T/\X\ for L^-^ = 0.30 and the experimental data of Ref.[7] 
(rescaled so that M^"^(0) = 0.35) using Jc/\X\ = 0.037 and 
the elastic constants given in the text. Inset is the pre- 
dicted average magnetization versus T/\X\ for a mixture of 
= 0.30 (5%) and L""/ = 0.33 (95%) components, with 
both magnetizations chosen to be positive just below Tc. 



moment drops as ^ jumps at the re-entrant JT tran- 
sition, Tj!^ is marked by a discontinuous change in 
M^^s(T). With decreasing L^^, both the orbital con- 
tribution Ml(T) and the magnitude of the change in 
M^^s(T) become smaller. For 1^ < L^^J < 0.324, the re- 
entrant JT transition causes the average magnetization 
to change sign, as shown in Fig.2(a) for 1%^ = 0.30. 

The quantitative agreement in Fig. 3 between the theo- 
retical prediction for L^^ =0.30 and the measurements of 
Ref.[7] for A=N(n-C4H9)4 is striking. Rather than pro- 
viding evidence for two compensation points, Ref.[7] indi- 
cates that a re-entrant JT transition occurs at Tj!^ ^ 42 
K. According to Fig. 2(b), a normal JT transition will be 

found at Tj!^ ^ 70K. Indeed, recent x-ray measurements 
[10] on the same compound confirm that the hexagonal 
symmetry present at room temperature is absent in the 
monoclinic lattice below 60 K. 

Measurements on several Fe(II)Fe(III) compounds [10, 
11] suggest that all MC compounds exhibit a small jump 
in the magnetization between Tc and Tcomp- The pre- 
dicted jumps in Fig. 2 (a) are much too large to explain 
those measurements. X-ray diffraction studies [12] re- 
veal that stacking faults in several compounds promote 
the coexistence of two phases: one with a six-layer repeat 
and the other with a two-layer repeat. The small mag- 
netization jumps observed [10, 11] in MC compounds are 
likely caused by a mixture of those two stacking types. 
As shown in the inset to Fig. 3, a mixture of two phases, 
type 1 with L^^ = 0.30 (5% of the sample) and type 



2 with L""/ = 0.33 (95% of the sample) [13], produces 
a small jump which is quite similar to the observations. 
The much larger jump observed by Tang et al. [7] may 
be caused by the greater fraction of type 1 {L^J ^ 0.30) 
stacking in their sample. 

While there is no JT distortion (^ = 0) when A < 0, 
^(T) remains nonzero down to T = when A > and 
L^^ < 0.260. So our model predicts that normal bimetal- 
lic oxalates with L^^ < 0.260 will not exhibit a disconti- 
nuity in the magnetization, which may explain why such 
a jump has never been observed in a normal compound 
[11]. Nevertheless, normal Fe(II)Fe(III) compounds with 
A > (so that the doublet remains lowest in energy) 

should manifest a normal JT transition at Tj!^ ^ 0.58|A| 
or about 85 K. 

But any non-Cs-symmetric cation like N(n-C4H9)4 will 
also induce a permanent distortion of the hexagonal lat- 
tice. Depending on the size and shape of the cation, 
this distortion can be local, weakly-correlated, or long- 
ranged. A non-Cs-symmetric potential can be included 
within our model by changing the off-diagonal terms in 
H_^^^ from ^ to ^ + ^0 • In the absence of spin-orbit cou- 
pling and a spontaneous JT distortion ^, the doublet 
splitting 6 is then 2^o- Since the anharmonic elastic po- 
tential proportional to favors ^ > 0, it acts to enhance 
the total distortion |C + Co| when > and to suppress 
the total distortion when < 0- For any nonzero ^o, 
there is a spontaneous JT distortion ^(T) 7^ at all tem- 
peratures due to the linear term of order — ^^o/|'^| in the 
free energy. 

As shown in Fig. 4 for L^^ = 0.3 and the elastic pa- 
rameters given earlier, increasing ^0 for a fixed L^^ de- 
creases while keeping Tcomp almost unchanged until 

Tjx < ^comp. When ^0 > 0.047|A| ^ 0.6 meV, the spon- 
taneous JT distortion ^ at T = jumps to a value near 
0.81 1 A| ~ 10 meV and the lower JT transition disap- 
pears. When ^0 < — 0.019|A| ^ -0.2 meV, the sponta- 
neous JT transition is eliminated but magnetic compen- 
sation survives until ^0 < — 0.43|A| ^ —5.4 meV, as seen 
in the inset to Fig. 4(a). Except in the narrow window 
> ^0 > — 0.019|A|, the permanent distortion ^0 and the 
electronic distortion ^ have the same sign. 

Clearly, the electronic JT distortion ^ will vanish if the 
Fe(II) singlet lies below the doublet with A < 0. So two 
conditions are required for C3 symmetry to remain un- 
broken at T = 0. First, the cation must preserve C3 sym- 
metry so that there is no permanent distortion ^o- Sec- 
ond, either A < or if A > 0, the average orbital angular 
momentum L^^ of the doublet must be sufficiently large 
to quench the spontaneous JT distortion ^ (L^-^ > 0.260 
in Fig. 2). Specific examples of symmetry breaking will 
be discussed in the next section. 
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FIG. 4: (Color online) The (a) average magnetization and (b) 
the JT mixing energy ^ normalized by |A| versus temperature 
T/|A| for several values of the permanent distortion ^o/|A| 
from -0.02 to 0.05 in increments of 0.01. Inset in (a) is the 
magnetization versus T/|A| for ^o/|A| between -0.5 and in 
increments of 0.1. These plots use L""/ = 0.3, Jc/\X\ = 0.037 
and the elastic constants given in the text. 



IV. FIRST-PRINCIPLES CALCULATIONS 

In order to estimate the orders of magnitude of the 
spontaneous JT distortion ^ and the permanent distor- 
tion we performed a series of first-principles calcula- 
tions within the framework of density-functional theory 
(DFT). We employed the local spin density approxima- 
tion (SDA) in the plane-wave-pseudopotential approach 
with the PBE [14] exchange correlation functional as 
implemented in the Quantum-ESPRESSO package [15]. 
We used Vanderbilt ultrasoft pseudopotentials [15-17] in- 
cluding, in the case of Fe [15, 17], d electrons in the va- 
lence and non-linear core corrections. An energy cut-off 
of 45 Ry was used. Since the orbitals at the Fermi level 



are very localized in this ionic solid, we used a single k- 
point. Spin orbit coupling was not included. In order 
to stabilize the electronic density, the calculations were 
performed with an electronic temperature of 0.02 Ry or 
270 meV. 

Predictions of the experimental electronic and mag- 
netic structure in highly localized systems are known 
to be significant challenges for most approximations of 
DFT. Indeed, we find that the energy difference between 
ferromagnetic and antiferromagnetic configurations is be- 
low the resolution of our theory. Moreover, the charge- 
density wave observed experimentally (which results in 
inter-penetrated Fe(II) and Fe(III) networks) is not found 
to be the ground state. This possibly signals the presence 
of strong self-interaction errors for the localized d orbitals 
[18]. The antiferromagnetic Fe(II) and Fe(III) ordering 
is stabilized only after enforcing the net spin within each 
unit cell to be equal to 1/2. 

Calculations were performed for several possible stack- 
ings (a6, aa6, abc) of the bimetallic Fe(II)-Fe(III) layers 
with A=N(n-C3H7)4 cations. The total energy differ- 
ences between these stackings were below the accuracy 
of the theory, consistent with the high degree of poly- 
morphism and the low energy cost for stacking faults 
suggested by the x-ray scattering experiments [12]. Since 
the relative energy difference A between the doublet and 
singlet is quite sensitive to the stacking of the bimetallic 
planes, the angular momentum L^-^ of the ground state 
may also depend on the stacking of the bimetallic planes. 

The doublet splittings and the order of magnitude of 
the JT distortions were estimated by focusing on a single 
bimetallic layer with different cations A. Within hexag- 
onal supercells, we studied the cations A=NH4 (1) and 
N(n-C3H7)4 (2). We also considered A=N(n-C3H7)4 (3) 
in a non-hexagonal supercell. 

(1) Although no oxalates contain the cation NH4 be- 
cause it is too small to stabilize the open honeycomb lat- 
tice, calculations on this system allow us to estimate the 
JT distortion in a Ca-symmetric environment. First, we 
relaxed the forces on all atoms constraining the symme- 
try to remain C3. We found that the Fermi level for the 
majority band lies at a doublet localized at the Fe(II) 
sites and occupied by a single electron. This configu- 
ration is consistent with our model, which requires the 
partial occupation of a doublet to explain the MC and 
the JT distortion. Second, we displaced the Fe(II) ion 
from the symmetric position on a grid of points r. From 
the ab-initio calculations, we extracted the total energy 
E{t) of the configuration and the energy splitting 5(r) 
induced in the partially-occupied doublet. Because the 
calculations where performed at an electronic tempera- 
ture (~ 270 meV) much larger than the splittings, both 
electronic levels of the doublet have an equal occupation 
of 1/2. At this large electronic temperature, no JT elec- 
tronic energy was gained and the symmetric point r = 
remained the position with minimum energy, as expected. 

We can, however, estimate the energy gain at T = 
as E{r) — 6{r)/2. This approximation assumes that the 
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total energy difference is given by the sum of the oc- 
cupied eigenvalues at zero temperature and neglects a 
small change in the electronic density [19]. With this ex- 
pression, we obtain the approximate Fe(II) position and 
doublet splitting at zero temperature. In absence of the 
spin-orbit coupling, we estimate that the Fe(II) moves 
0.03 A (the nearest-neighbor Fe(II)-Fe(III) distance is 
about 5.4 A) with an energy gain of 2 meV and a doublet 
spitting of 8 meV. This JT splitting is quite close to the 
spin-orbit coupling {XlL^-^ S ~ 7.5 meV with L^-^ = 0.3. 

(2) Moderately-sized cations such as N(n-C3H7)4 are 
themselves non-Ca-symmetric and break the C3 symme- 
try of the crystal-field potential at the Fe(II) sites. While 
it is possible to construct a N(n- 03117)4 isomer with three 
C3H7 radicals below the oxalate plane, thereby preserv- 
ing C3 symmetry, the remaining propyl chain (oriented 
towards the hexagonal hole in the oxalates plane) must 
break C3 symmetry. This symmetry violation produces a 
crystal-field splitting 2^o of the doublet and a permanent 
distortion of the open honeycomb lattice. 

Calculations were performed in a hexagonal unit-cell 
containing a single A=N(n-C3H7)4 cation with periodic 
boundary conditions, corresponding to to an ordered con- 
figuration where every N(n-C3H7)4 cation is oriented in 
the same way. After relaxing the positions of the atoms, 
we obtained an intrinsic doublet splitting 2^o of about 10 
meV at an electronic temperature of 270 meV. 

These calculations also indicate that the intrinsic dis- 
tortion introduced by the cation will be increased at low 
temperature by the electronic energy gain of the Fe(II) 
JT distortion. So in the absence of spin-orbit coupling, 
the parameters ^ and ^0 of our model have the same sign 
for this cation. 

(3) The molecule N(n-C3H7)4 is just small enough 
to allow individual cations to rotate independently of 
each other. Hence, a uniform distortion of the crystal 
is possible. Such a uniform distortion of the open honey- 
comb structure has been observed in bimetallic oxalates 
with radical cations [20]. For an ordered arrangement 
of A=N(n- 03117)4 cations and allowing the supercell to 
break hexagonal symmetry, we obtain a doublet splitting 
^ ~ 20 meV, which should increase to about 30 meV 
when the electronic temperature approaches zero. 

By contrast, larger cations such as A=N(n-C4H9)4 
studied in Ref.[7] are unable to independently rotate 
within each unit cell. Because they are locked into 
a highly-disordered configuration during synthesis, such 
cations will only break the local C3 symmetry around 
each Fe(II) ion. Since the spin-orbit coupling competes 
with the small, local lattice distortions, magnetic com- 
pensation is still possible when the cations are large 



and disordered. This leads to an interesting conjecture: 
due to their ability to spatially order, smaller cations 
like N(n- 03117)4 may be more eff'ective at enhancing the 
departure from O3 symmetry and suppressing magnetic 
compensation than larger cations like N(n-04H9)4. 
V. DISCUSSION 

Our first-principles calculations imply that large 
cations like N(n-04H9)4 will induce a substantial, local 
distortion ^0 of the hexagonal lattice. For smaller cations 
like N(n-03H7)4, this distortion may propagate through- 
out the bimetallic layer because the cations are able to 
reach a global minimum of the free energy, which is an 
ordered state at low temperatures. For an ordered con- 
figuration of N(n-03H7)4 cations with |^o| ^5 meV, the 
value suggested by our first-principles results. Fig. 4 im- 
plies that the re-entrant JT transition will be absent and 
any magnetic compensation will be very weak. Indeed, 
N(n-03H7)4[Fe(II)Fe(III)ox3] bimetallic oxalates [5] show 
no signs of a JT transition or magnetic compensation be- 
low Tc. 

With the spin-orbit coupling set to zero (L^-^ = 0), the 
doublet splitting S = 2^(0) ^ 21 meV obtained using 
the parameters of Fig. 3 is more than twice as large as 
that obtained from first-principles calculations in a O3- 
symmetric environment (see 1 above). So it is likely that 
the organic cation N(n-04H9)4 plays a significant role in 
breaking the local O3 symmetry and enhancing the dou- 
blet splitting. A more sophisticated description of the 
experimental measurements may be possible once addi- 
tional information about the atomic structure becomes 
available. 

To summarize, we have provided strong evidence 
for the existence of a re-entrant JT transition in the 
Fe(II)Fe(III) bimetallic oxalates. Observation of the JT 

distortion between t]^ ^ 42 K and Tj^^ ^ 70 K would 
provide unambiguous support for the predicted multiplet 
splitting in MO compounds. We hope that this work will 
inspire systematic x-ray scattering measurements that 
will verify the predictions made in this paper, includ- 
ing the long-range ordering of small, non-03-symmetric 
cations in this important class of layered, molecule-based 
magnets. 
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